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In a recent paper we presented a calculation of NNLL virtual corrections to the forward- 
backward asymmetries in 6 — Xsi~^i~ decay. That result does not include bremsstrahlung 
corrections which are free from infrared and collinear singularities. In the present paper we 
include the remaining 0{as) bremsstrahlung corrections to the forward-backward asym- 
metries in 6 — i> Xsi~^i~ decay. The numerical effect of the calculated contributions is 
found to be below 1%. 



1 Introduction 



Rare B-decays are known to provide a unique source of information about the physics 
at the scales of several hundred GeV. In the standard model (SM) all these decays proceed 
through loop diagrams and are suppressed. Thus in the extensions of the SM the con- 
tributions from the 'new' sources of flavor-violation can be comparable to or even larger 
than the SM contribution. Therefore experimental information on rare decays can be 
used to test the SM at the one-loop level or to put constraints on its extensions. 

During the last decade the experimental and theoretical efforts were concentrated on 
the 6 — s> S7 mediated decays. Available experimental data already provides stringent 
constraints on certain extensions of the SM. In 2002, also the exclusive B Kgix^ n~ 
decay mode was measured by the BELLE collaboration [1]. This measurement was con- 
firmed soon by the BABAR Collaboration [2]. Recently, also the first measurement of the 
branching ratio in the inclusive decay B Xsi^i^ has been reported by the BELLE Col- 
laboration [3]. The results of these first mcasurmcnts are compatible with SM predictions 
though more statistics will be needed for more decisive conclusions. It is expected that 
precise measurement of kinematical distributions for the B — > Xsi'^i~ decay combined 
with data on S — > Xg^ will significantly tighten the constraints on the extensions of the 
standard model [4]. 

^From the theoretical point of view the description of B ^ X^i^i^ decay is prob- 
lematic because of the long-distance contributions from intermediate cc resonant states. 
When the invariant mass y/s of lepton pair is close to the mass of resonance, only model 
dependent predictions for these long distance contributions are available today. However 
for the region 0.05 < s = s/ml < 0.25 the nonperturbative effects are estimated to be be- 
low 10% [5]- [10] and the differential decay rate for B — >■ Xs£~^£~ can be well approximated 
by HQET corrected short distance contribution. 

The next-to-leading logarithmic (NLL) calculation for B — > Xs£'^i~ has been per- 
formed quite long ago in [11, 12]. However those results are known to suffer from a 
relatively large (±16%) dependence on the matching scale fiw- The NNLL corrections to 
the Wilson coefficients eliminate the matching scale dependence to a large extent [13], but 
leave a ±13%-dependence on the renormalization scale /ib, which is of order 0(mb). To 
further improve the theoretical prediction, 0(0;^) virtual and bremsstrahlung corrections 
have been calculated [14]- [16]. As a result the renormahzation scale dependence reduced 
by a factor of 2. 

It is well known that the measurement of the forward-backward asymmetries for de- 
cay B — > Xgi^i^ can be used, in combination with the measurement of B ^ ^sl, to 
perform a so-called model independent test of the standard model [4], [17], [18]. The dou- 
ble differential decay width dr/[dsdcos9] and the forward-backward asymmetries have 
recently been calculated with NNLL precision [19], [20]. Those calculations included one- 
loop virtual 0{as) corrections associated with the operators O7, Og and Oiq as well as 
the corresponding bremsstrahlung corrections that are necessary to cancel the infrared 
and mass singularities. It was found that NNLL corrections drastically reduce the renor- 
malization scale dependence of forward-backward asymmetries. In the present paper we 
complete the calculation of the NNLL calculation for the forward-backward asymmetries, 
presenting the full results for the bremsstrahlung corrections associated with the operators 
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Oi, O2, Os which were omitted in [19], [20]. 

The paper is organized as follows. In section 2 we briefly describe the theoretical 
framework. In section 3 the analytical results for the forward-backward asymmetries in 
b — > Xs£'^£~ decay are presented. In section 4 we discuss the technical details of the 
calculations and give phenomenological analysis for the forward-backward asymmetries 
in 6 — >• Xgi'^£~ decay . 



2 The Theoretical Framework 

The most efficient tool for studies on weak B meson decays is the effective Hamiltonian 
technique. For the specific channels b sl'^l" {I — ii,e), the effective Hamiltonian is of 
the form 

AC 

neff ^ -^v*Vt,j:QO, (1) 

v2 



where we have omitted the contribution proportional to the small CKM factor V*gVub- 
The dimension six effective operators can be chosen as 

Oi = {sl1,T-Cl){clYT%l), O2 = (sl7mCl)(cl7^&l), 

Q q 

9 9 

O7 = \mb{sLa''''bR)F^,, Os = -mb{sLa^''bR)F^,, (2) 

9s 9s 
2 2 

09 = ^ (sl7m&l) iV) , Oio = ^ {sLlM (^7m75^) 

9s 9s 

where the subscripts L and R refer to left- and right- handed components of the fermion 
fields. In the following it is convenient to use the related operators O7, ...,Oio, defined 
according to 

Oj^^O, (j = 7,...,10) (3) 
with the corresponding Wilson coefficients 

Cj^^Cj (j = 7,...,10). (4) 

We refer to [13] and [15] for details. 



3 NNLL results for the forward- backward asymme- 
tries 

We start introducing the forward-backward asymmetries. Often the so-called normal- 
ized and un-normalized forward-backward asymmetries are considered which are defined 
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by 

and 

Jd^r(& ^ Xd^nidsdz\sgn{z)dz ^ ^ 
^^^^'^ = r(^ X.ePe) 

respectively. 

The double differential decay witdh can be written in the form [19]: 

d'r{b ^ x^in-) _ fae.m.y Glml^,i,\V*Vtbl _ 2 



dsdz V 47r / 4871^ 



+?[(! + + ^(1 - ^^)] la^^^r (l + — /77(^, ^)) - 3^^i?e((79^^^Cri(^*) (l + ^/9io(^) 

+6i?e(C',^^^C'9^^^*) (l + ^/79(s, ^)) - 6^i?e(C7^^^C,t^*) (l + ^/7io(5)) } (7) 



where z — cos{9) and the effective Wilson coefficients C^^-^ , Cg-^^ and C^q-^ are given by 

Cg^^^ = ^ + Tgh{ml, s) + Ugh{l, s) + Wgh{0, s) 

«.(/^) p(o)^(9)(^) ^ (5) + (5)], (8) 



<-^io — ^10, 

with ml = ml/ml The quantities C}°\ \ Aj, Ag, Aio, Tg, Ug and Wg are Wilson 
coefficients or linear combinations thereof. Their analytic expressions can be found in [14]- 
[19]. In Eq. (7) the term ^^^'^dsdz"^^^ finite bremsstrahlung corrections 

that have not been considered in [19]. In this paper our goal will be to investigate the 
impact of that term on forward-backward asymmetries. The general distribution (7) will 
be investigated elsewhere [21]. 

In the numerator, both asymmetries involve the same integral that can be expressed 
as follows: 



7-1 dsdz ^9n{z}az -y^^) 4g^3 



2 



X 



-3sRe (Cl^Ctt) (l + ^/9io(^)) - QRe {Cr"Ctt) (l + '^fno{S] 



(9) 
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^810 — 



rrerence or the matrix elements or 
the present paper. Their exphcit expressions read: 

— ((3((1 - Vi)'(23 - eVi - s) + 4(1 -s){7 + s) ln(l + VI) 
d(1 — s)"^ 

+ 2s(l + s- ln(s)) ln(s))) + 2(-37r2(l + 2s) + 6(3 - s)sln(2 - VI) 
- 36{l + 2s)Li2[-VI] - W^^(6(2(-3 + s)sarctan 



to 

,08 



2 + VI' 

[V^ 



+ 27rln(2-Vi) — arctan 

'VIV4: - S 



((-3 + s)s + 41n(2- Vi)) 



— arctan 



2-s 



- 2Re{iLi2[^V4^{l - s)VI + ^^—^ 



(( Q^-^- 1 (-w^A-D CT- A-'^ + ^V^- s + Vsh)VI.. 
((-3 + s)s - ln(s)) + ARe[iLi2[- J) 

iv4 — s — vs 



■]))))) 
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C210 



+ 



where 



{(4(l-l)(l + ^)- 



2J{s - 'u;2)2(^/;(3 + w) - s(l - w)) 



Is (s — w)(l — s)^ 

+ {2 + 5w + 2w^ + s{3 + Aw)) In 



s — w 



2\2 



2w 



[s — w) 



S^(l +U))2 -4s 



X (w(2 - w) - s(6 - 5w))(ln(l + w - s(3 - w) + (1 - s)^(l + w)^ - 4s) 
- ln(s(l - 3w) + w^{l + w) + ^(s- ^2)2^(1 + w)^ - 4s)))A2 



^23 



- (2(l-s)(l + 2w) - 



2 J(s - w^)'^{w{2 + w)- s(l - w)) 



+ 2(s(l + 2w) +w(2 + w))ln 
4(1 - w){s - w) 



+ + ^(s - w'^y 

2w 



, (ln(l + w - s(3 - w) + (1 - s)V(l + w)2 - 4s) 

^(1 + ^)2-45 

ln(s(l - 3w) + w2(l + w) + J{s - w2)2 J(l + w)2 - As)))A2'j}dw 



A^,3 = -2 + 



Ai,, = 2 



w — s 

s 



The functions G-i(t) and G'o(t) can be found in [16] (Eqn. (30), (31)). 
The functions ^210 and t^io are plotted in Fig.l. 



(11) 



(12) 
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4 The details of the calculation and numerical results 



As in [19] we will work in the rest frame of the lepton pair. The corresponding formulae 
have been derived in [19]. For the 3- momenta of b-quark (pb), {h) and gluon (r) in 
that frame of reference we have: 

Pb = (Ip&I, 0, 0), k = {^cose, ^sine, 0), 
f = {ErCosOi, ErSin9iCos92, ErSin9isin92) , 

where 

2Eb^/s -2Er^/s + 2ErEb - s-ml + w? 



cos9i — 



s 



2Er\pb\ 



Eb is the energy of the b quark and Er is the energy of gluon. Then the formula for double 
differential decay width reduces to 

d^T{b sgi+r) rUbS 



dsdz (27r)62^ 



X 



j |M|2(1 - Z2^)~^/^dz2dErdEb (13) 



where |Mp is the squared matrix element, summed and averaged over spins and colors 
of the particles in the final and initial states, respectively, s — s/ml and z — cos9, 
Z2 = cos92- The boundaries for the integration variables are 

pmm ^ rub'^ + s- 2Eb^s - rris^ < ^ < m^^ + g - 2Eb^s - ^ 
2{Eb^\Vb\-^s) - 2{Eb-\vb\-^s) 

mb<Eb< " ^ ^ -\<Z2<\ (14) 

l\J s 

As we have mentioned above the forward-backward asymmetries arise only from inter- 
ference of Oio matrix element with Oi,i = 1, ..,9^. The infrared infinite bremsstrahlung 
corrections coming from interference of Oiq with O7 and Og (along with the correspond- 
ing virtual corrections) have been taken into account by introducing functions /710 and 
/910 which have been calculated in our previous paper. The remaining bremsstrahlung 
corrections are infrared safe and are coming from the interference of Og and Oi, O2 with 
Oio (see Fig. 2 and Fig.3 for the contributing Feinman diagrams). 

The calculation of and Oio interference term is relatively easy and can be performed 
analytically. The calculation for Oi,2 and Oio interference terms is more complicated and 
will be described in some dctailcs below. 

The bremsstrahlung corrections involve the matrix elements associated with the two 
diagrams in Fig. 3. Their sum, J^^, is given by 



— G 9s Qu 

16 TT^ 



E{a, P, r) Ah + E{a, /3, q) Aig - E{P, r, q)—Ki2^ 

qr 



E{a, r, q)^ Aise - E{I3, r, q)^-^ M27 
qr qr 



(15) 



^As in our previous papers we will systematically ignore the 0{as) contributions coming from matrix 
elements of operators Oj, i = 3, 4, 5, 6 as their Wilson coefficients are very small. 
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where q and r denote the momenta of the virtual photon and of the gluon, respectively. 
The index a will be contracted with the photon propagator, whereas will contracted 
with the polarization vector e^(r) of the gluon [15]. The matrix E{a,P,r) is defined as 

E{a, f3, r) = ^(7a7/3^ - hpla)- (16) 

Due to Ward identities, the quantities Aik are not independent of one another. Namely, 
we have 

q^ - 

Ais = Ai23 + — Ai27 ; Aie = Aise • (17) 
qr 

As in addition Az26 = — Ai23, the bremsstrahlung matrix elements depend on Ai23 and 
Ai27, only. In 0? = 4 dimensions we find 



Az23 = 8 (qr) ['dx dy ""^^^^ , At27 = 8 (qr) ['dx dy yJ^—lt , (18) 
Jo L> Jo C 



^0 

where 

C^ml-2xy{l- y){qr) - q^ y {1 - y) - iS. (19) 

It is easy to notice that dependence of the matrix element M on Z2 is polynomial so the 
integral over Z2 in Eq.(13) is straightforward. To deal with the remaining two dimensional 

integrals over Er = Er/mb and Ei, = Eb/rrib it is useful to introduce a new integration 
variable w instead of Ej. defined by w = 2\f§Ej. + s .The integration limits are then given 

by 

s^-w^ 1 + s 



and w e [s, 1]. (20) 



For a fixed value of s, the quantities Ai23 a-nd Ai27 depend only on the scalar product qr, 
which is given by {w — s) ml/2. The integration over Ef, then turns out to be of rational 
kind and can be performed analytically. The remaining integration over w, however, is 
more complicated and is done numerically. 

We now proceed to the investigation of the numerical impact of the finite bremsstrahlung 
corrections on the forward-backward asymmetries. In Fig. 4 we show the contribution of 
the finite bremsstrahlung corrections to the unnormalized forward-backward asymmetry 
AAfb{s)- Dashed-dot lines show the contribution from interference of matrix elements 
of the operators Og and Oiq (fi =2.5 GeV (uppermost curve), 5 GeV (middle curve), 
and 10 GeV (lower curve)), dashed lines show the contribution from interference of ma- 
trix elements of the operators Oi,2 and Oio {ji =2.5 GeV (lower curve), 5 GeV (middle 
curve), and 10 GeV (uppermost curve) and ), solid lines show the sum of previous two 
terms (/x =2.5 GeV (lower curve), 5 GeV (middle curve), and 10 GeV (uppermost curve)), 
mc/mh=0.29. Wc see that contributions coming from Oi,2 — Oiq and Os — Oio interfer- 
ence terms partially cancel each other. For that reason the contribution of new terms 
to the forward backward asymmetries is numerically small. In Fig. 5 we combine the 
new corrections with the previous results for unnormalized forward-backward asymmetry. 
The solid lines show Apb{s) including new corrections for three values of renormalization 
scale {/I— 2.5,5 and 10 GeV) and rfic = 0.29. The dashed lines represent the corresponding 
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results without finite bremsstrahlung corrections.. In Fig. 6 we give the correspond- 
ing curves for normalized forward-backward asymmetry Afb{s). The new contribution 
is numerically small: at s = we find Ap^^^{0) — 2.304 ± 0.109 while without finite 
bremsstrahlung corrections we have Ap^^^{0) = 2.292 ±0.111. For the position sq where 
forward-backward asymmetry is zero we find s^^^^ = 0.1620±0.0016 while without finite 
bremsstrahlung corrections we have Sq^^^ = 0.1615 ± 0.0015. 

It is necessary to mention that for some extensions of the standard model (see for 
instance, [22]) with opposite sign (and larger values) of C^-^^ the contribution of the finite 
bremsstrahlung corrections to the forward-backward asymmetries can be sizable. 

To conclude, we have calculated the finite bremsstrahlung corrections to the forward- 
backward asymmetries for b — > Xsi~^i~ decay. We found that the numerical impact of 
the new corrections on the forward-backward asymmetries in the standard model is less 
than 1%, but for some extensions of the standard model it can be larger. 



ACKNOWLEDGMENTS 

We thank C. Greub for useful discussions. The work was partially supported by 
NFSAT-PH 095-02 (CRDF 12050) and SCOPES programs. 



Table 1: Coefficients appearing Eq. (8) for fj, — 2.5 GeV, ji — h GeV and = 10 GeV. 

For as{ij) (in the MS scheme) we used the two-loop expression with five fiavors and 
(^si'iTT'z) = 0.119. The entries correspond to the pole top quark mass rrit — 174 GeV. The 
superscript (0) refers to lowest order quantities. 







2.5 GeV 


5 GeV 


10 GeV 


Cts 




0.2675 


0.2150 


0.1800 






-0.6960 


-0.4860 


-0.3257 






1.0458 


1.0235 


1.0109 






(-0.3605, 0.0322) 


(-0.3207, 0.0185) 


(-0.2873, 0.0089) 


A^ 




-0.1688 


-0.1535 


-0.1403 






(4.2407, - 0.1700) 


(4.1283, 0.0130) 


(4.1311, 0.1549) 


(Tf, 




(0.1148, 0.2791) 


(0.3744, 0.2512) 


(0.5763, 0.2307) 


(f/f, 




(0.0455, 0.0228) 


(0.0326, 0.0153) 


(0.0224, 0.0105) 






(0.0440, 0.0162) 


(0.0321, 0.0117) 


(0.0223, 0.0085) 


(4?, 


AS?) 


(-4.3731, 0.1349) 


(-4.3731, 0.1349) 


(-4.3731, 0.1349) 
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Figure 1: Functions tgio (upper curve) and ^210 
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Figure 2: Bremsstrahlung diagrams induced by Oi and O2. 




Figure 3: Bremsstrahlung diagrams induced by O7 and Og, Og, Oio 
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Figure 4: The contribution of finite bremsstrafilung corrections to tfie unnormalized 
forward-backward asymmetry. Dashed-dot lines show the contribution from interference 
of matrix elements of operators Og and Oiq {fi =2.5 GeV (uppermost curve), 5 GeV 
(middle curve), and 10 GeV (lower curve)), dashed lines show the contribution from in- 
terference of matrix elements of operators (9i_2 and Oio {fj, —2.5 GeV (lower curve), 5 GeV 
(middle curve), and 10 GeV (uppermost curve) and ), solid lines show the sum of finite 
bremsstrahlung contributions (yU =2.5 GeV (lower curve), 5 GeV (middle curve), and 10 
GeV (uppermost curve)), rric/mb— 0.29. 
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Figure 5: The solid curves show s and ii (for /x=2.5, 5 and 10 GeV) dependence 
of unnormahzed forward-backward asymmetry in NNLL approximation including finite 
bremsstrahlung corrections while dashed lines show the corresponding results without new 
corrections. 
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Figure 6: The solid curves show s and n (for /x=2.5, 5 and 10 GeV) dependence 
of normahzed forward-backward asymmetry in NNLL approximation including finite 
bremsstrahlung corrections while dashed lines show the corresponding results without 
new corrections. 
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